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ABSTRACT 

We revisit a method to obtain upper limits on the jet matter content combining synchrotron 
self-absorption constraints and the large scale bubble energy. We use both X-ray observations, 
which give limits on the jet power from the energies and timescales of bubbles found in 
clusters of galaxies, and radio observations, which give limits on the magnetic field in the 
jets. Combining the two imposes constraints on the particle number density, and hence the 
jet content. Out of a sample of clusters which have clear radio bubbles, there are only two 
which have sufficient resolution in the radio images to give significant constraints, under the 
assumption that the jets are fairly steady. The results for M87 and Perseus indicate that the 
radio emitting region of the jet is electron-positron dominated, assuming that the minimum of 
the electron energy distribution, 7min ~ !■ 
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1 INTRODUCTION 

The radio synchrotron emission from jets associated with active 
galactic nuclei (AGN) indicate that they must contain relativistic 
leptons and magnetic fields. The matter content of these jets, how- 
ever, is still a mystery (e.g. Homan 2005). Assuming that these 
jets are electrically neutral, are they electron-positron ("light") or 
electron-proton ("heavy")? This paper aims to determine whether 
the positively charged particles are positrons or protons, in the as- 
sumption that the magnetic field is not energetically dominant. 

Previous attempts to determine the matter content of the jets 
have used Synchrotron Self-Absorption (SSA) arguments (e.g. 
Celotti & Fabian 1993; Reynolds et al. 1996). Measurements of the 
core radio flux and the size of the emitting region lead to constraints 
on the number density of the relativistic particles responsible for 
the emission, n, and on the magnetic field present in the jet, B. The 
effect of any protons in the jet cannot be determined from radia- 
tive information. As a result the large scale properties of the source 
are required to place independent estimates on the total power. In 
the past this has been difficult because the total energy in the radio 
lobes has not been known; equipartition between the energy in the 
magnetic field and the radio-emitting particles has been assumed in 
order to be able to estimate the energy present (Burbidge 1959). 

With the advent of the Chandra and XMM-Newton X-ray ob- 
servatories the interaction of radio sources in the centres of clusters 
of galaxies with the surrounding Intra-Cluster Medium (ICM) has 
been observed in great detail. One of the first observed was the 
Perseus Cluster, where decrements in the X-ray emission from the 
thermal ICM were first seen by ROSAT (Bohringer et al. 1993). 
These have been interpreted as "bubbles" blown in the ICM by the 
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central radio source. Strong support for this interpretation came 
with Chandra imaging and spectral information (Fabian et al. 
2000). Subsequently many other examples of such bubbles have 
been found (e.g. Hydra A, McNamara et al. 2000; A2052, Blanton 
et al. 2001; A2199, Johnstone et al. 2002; Centaurus, Sanders & 
Fabian 2002). 

As these radio lobes are embedded in the ICM, their en- 
ergy content can be more accurately determined using arguments 
of pressure balance between the thermal ICM and the relativistic 
plasma within the bubbles, as well as bubble dynamics. As a re- 
sult, the connection between the kinetic luminosity of the jet and 
the energy in the bubble can be calculated. The potential of this ap- 
proach is that it does not require assumptions on equipartition for 
the estimate of the jet power. 

In this work we combine the information from X-ray and 
VLBI observations of radio sources embedded in clusters of galax- 
ies, where clear bubbles are visible in the X-ray images. We apply 
the calculations performed for the jet in M87 by Reynolds et al. 
(1996) to calculate the matter content of the jet present in this and 
other sources. In essence the method determines the lepton content 
of the base of the jet from synchrotron emission and absorption, 
and the total particle energy content from power required to cre- 
ate the observed bubbles. The comparison between the two, sets 
constraints on the particle content. 

In Section 2 we outline the models and assumptions used in 
the calculations. The source parameters and the results are pre- 
sented in Sections 3 and 4. The assumed lepton energy distribu- 
tion is discussed in Section 5, and further complications of this 
method and comparisons with other methods are discussed in Sec- 
tion 6. The implications of this work on the particle content of radio 
lobes is investigated in Section 7. We use Ho = TOkms" 1 Mpc" 1 
throughout. 
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2 THE MODEL 

This section briefly summarises the calculations presented in 
Reynolds et al. (1996) (following Marscher 1987) in order to 
demonstrate the method and define the parameters. 

2.1 Limits on the Magnetic Field from Synchrotron 
Self-Absorption 

With the assumption that the relativistic leptons have a distribution 
of Lorentz factors N(-y) = iVo7~ p = N j~ (2a+1 \ between 7 min 
and 7 m ax, the spectrum in the optically thin regime is S v oc v~ a . 
The corresponding relativistic lepton number density, n, is given 

by 



Hnui 

n = N(j)dj 
2an = -No [l~ 2a ] 



2a] 7max 
7min 

Therefore, for 7 max >> 7 m i n and a > 0.5, the lepton number 
density is 

No 2a m 

Initially we assume 7 m in = 1 for the following calculations. For a 
discussion of the effect of varying 7 m in see Section 5. 

The synchrotron flux in the optically thick (self-absorbed) re- 
gion is independent of the particle density, resulting in an estimate 
on the magnetic field of 



b < irr 8 b(a)e^s~ 2 



l + z 



G. 



(2) 



where the u m and S m are the frequency (in GHz) and flux den- 
sity (in Jy) at spectral turnover. 9 d is the angular diameter of the 
source in milliarcseconds (mas). As the angular diameter may be 
an upper limit (the core may not be fully resolved) this causes the 
estimate on B to be an upper limit. In some cases the source is el- 
liptical, and in these cases the average 9d = VQJK, where # a and 
9b are the corresponding angular diameters of the ellipse. 8 is the 
relativistic Doppler factor, defined as 8 = l/r(l — /3 cos <f>) where 
r and [3 are the Lorentz factor and v/c for the bulk motion of the 
jet respectively and <f> is the angle between the line of sight and the 
jet axis. The function, 6(a), is tabulated in Marscher (1987) and we 
have interpolated between the values given where appropriate. 

The observation of compact radio sources with flat spectra, 
Qobs ~ 0, has been interpreted as the superposition of different 
SSA components each peaking at different frequencies (Blandford 
& Konigl 1979). Hence the observations at a given v basically mea- 
sure the flux density and the size of the component of the jet which 
is becoming self absorbed at u (in the observer's frame). 

2.2 Limits on the Magnetic Field and Number Density 

The jet becomes self-absorbed at an observed frequency, z/ m , where 
T sy n(Vm,r) = ft(f m )X = 1; where n(y) is the synchrotron ab- 
sorption coefficient, X is the path length of the line of sight through 
the jet and r the jet cross-section radius. From Reynolds et al. 
(1996) and references therein 

= 3 (Q+1) v ^ : g(a)e 2 7Yo ^ ( 3 /2+a) ^ (5/2+a))5(5/ 2 +a )^ 
m 8m c c B m 

(3) 

where u B is the cyclotron frequency and g(a) is the product of 
gamma functions (of order unity for the considered range of a). 



This expression for n(y) is valid for v 3> i/ m i n ~ 7min^B, 
where Vmin is the low-energy cut-off in the spectrum correspond- 
ing to the low-energy cut-off in the lepton energy distribution. In 
this case the self-absorption only depends on the normalisation of 
the relativistic lepton distribution and the magnetic field. 

The path length is X = 2r /8, using the relativistic transfor- 
mations for a cylindrical geometry. Combining this with Equations 
1 and 3, places a lower limit in the B — n plane for radiation at 
frequency i/ m to be self-absorbed in the source, 
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2.3 Kinetic Luminosity 



The kinetic luminosity of the jet depends on the type of parti- 
cles present, as well as their energy and number density. In the 
assumption that all of the energy contained within the jet re- 
sults in the creation and expansion of the radio bubbles observed 
within clusters, then the power required to create these bubbles 
is an estimator of the (average) kinetic luminosity of the jet, 
i.e. Lk = -Ebubbie/tbubbio, where ^bubble is the creation time of 
the bubble. The simplest estimate on the energy required is that 
^bubble = pV, where V is the volume of the bubble, and p is the 
pressure of the surrounding intra-cluster gas. Taking into account 
any internal energy of the bubble results in 

7R 



E h 



r pV, 



7R _ 1 

where 7r is the ratio of specific heat capacities, which for a rela- 
tivistic gas is 4/3, resulting in Sbubbic = 4pV. Whether the energy 
contained within the bubble is pV, 4pV or some other multiple 
of pV is currently uncertain. For example, investigating the weak 
shock surrounding the bubbles in the Perseus cluster, Fabian et al. 
(2006) find that the energy of the post-shock gas is around 2pV. 
We use iJbubbic = pV, although using 4pV does not change the 
results significantly. 

In our analysis we assume that the radio bubbles are in pres- 
sure equilibrium with their surroundings, whereas Reynolds et al. 
(1996) assumed that they were over-pressured by a factor of ~ 3. 
Using their source parameters we recover the limits they place on 
the matter content of the M87 jet (see their Fig. 1). 

There are a number of estimates on the bubble timescales. The 
most appropriate one for these young bubbles which are (presum- 
ably) still being inflated by their jet is the sound speed timescale, 
^bubble = 2i?bubbio/c s , where ^bubble is the bubble radius and c s 
is the local sound speed, following Dunn & Fabian (2004). There 
are no indications for strong shocks surrounding the bubbles, which 
implies that the bubble edges are travelling at less than the local 
sound speed. It is possible, however, that the bubbles do not grow 
smoothly, but in fits and starts (Fabian et al. 2005), and as such this 
timescale is not a good estimate for the age of the bubble. For fur- 
ther discussion on timescales relevant to the evolution of bubbles in 
clusters see Churazov et al. (2000); Dunn & Fabian (2004); Dunn 
et al. (2005). 

As mentioned, we assume that the jet is particle dominated 
(i.e. we neglect the magnetic field contribution) and for simplicity, 
following other work (e.g. Sikora & Madejski 2000), assume that 
the protons are "cold". Thus the jet kinetic luminosity, including 
the advected energy, is 



L K « r 2 (3TTr(Z) 2 nm e c 3 



^«7>-l) + L r ^(l + fc a ) 



(5) 
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Table 1. : Cluster Values and Notes 



Bubble 




-Ebubblc 


^bubble 




Degrees 


10 57 erg 


10 7 yr 




M87, 2 = 


= 0.004 




Jet 


15.0 ±5.0 


0.28 ±0.14 


0.59 ± 0.03 


Counter- Jet 


i £ n _L c n 
15. ± 5.0 


f\ qo _l_ n on 
O.oy ± 0.20 


o A o _l_ f\ no 
0.48 ± 0.02 


Average 




0.33 ±0.12 


0.53 ± 0.02 


L K = 


2.0 X 10 42 ergs" 


■ 1 






Perseus, z 


= 0.018 




Northern-Inner 


40.0 ± 10.0 


4.97 ±2.20 


0.98 ±0.18 


Southern-Inner 


40.0 ± 10.0 


4.32 ±1.98 


0.93 ±0.18 


Average 




4.64 ± 1.48 


0.95 ±0.13 


L K = 


1.5 X 10 43 ergs" 


■ l 





REFERENCES: The energies and timescales for M87 and Perseus are cal- 
culated from the bubble radii presented in Allen et al. (2006) and Dunn 
et al. (2006) respectively. The angles to the line of sight are from Biretta 
et al. (1999) and Walker et al. (1994). Reynolds et al. (1996) used L K = 
lO^ergs- 1 . 



where fc a takes into account the effect of hadrons on the rest-mass 
energy (adapted from Schwartz et al. 2006). For electron-positron 
jets fc a = 0. In electron-proton jets we assume that there is one pro- 
ton per electron, so fc a = m p /m e . The internal energy of the lep- 
tons is expressed as ((7) — l)nm c c 2 ~ (^)nm e c 2 for (7) >> 1. 
If a = 0.5 then the average particle Lorentz factor is 



(7) = Tminln 




For a more general case, where 7 max >> 7 m i n , 

, \ 2a 7 gg n ( 1 1 \ „, 2a 

{V ~ 1 - 2a { Jf-i) (2a-i) I ~ 2a - l 7 ™"' 

\ ' max 7min / 

Using the bubble power and Equation 5 estimates on the num- 
ber density can be calculated for the two different jet compositions. 
These expected number densities are compared to the limits im- 
posed on the B — n plane by Equations 2 and 4. 



3 SOURCE PARAMETERS 

Here, as a result of the range of spectral indices measured, we as- 
sume that each SSA component has the same particle energy dis- 
tribution, with p = 2.4 and so the same (optically thin) spectral 
index a = 0.7, as this falls in the middle of the range expected for 
optically thin spectra (0.5 < a < 1.0). 

The core fluxes presented in this work are those which place 
some useful limit on the matter content of the jets. There are many 
other, mainly older, measurements which result in lower limits on 
the number densities inferred which are too low to allow any con- 
straints to be placed on the matter content. 

3.1 M87 

We include the VLBI and Lk parameters adopted by Reynolds 
et al. (1996) (using a = 0.5) to allow a comparison to be made 
between their values and the values based on updated observations. 
The jet parameters used Reynolds et al. (1996) were F > 3 and 



3 

4> ~ 30°, based on the observations in Biretta et al. (1991) and 
Biretta et al. (1995). 

The apparent jet speeds measured by Biretta et al. (1995) 
ranged from 2c close to the core to ~ 0.5c for more distant com- 
ponents. Recent observations report a range of apparent speeds; 
Kellermann et al. (2004) measure 0.04 ± 0.02c for a compo- 
nent 6 mas from the core; Ly et al. (2004) place lower limits of 
0.25 - 0.40c on the jet motion and Dodson et al. (2006) find lit- 
tle evidence for any motion of the components close (< 150 mas) 
to the core. A series of HST images resulted in measurements of 
apparent speeds of 2.6 — 6c for components between 0.87" and 
6" from the core (Biretta et al. 1999), with the innermost resolved 
component at 160 mas from the core having an apparent speed of 
0.63 ± 0.27c. 

Clearly the superluminal velocities measured might not be as- 
sociated with the bulk flow of the jet, but correspond to pattern 
speeds (e.g. shocks). However, if the observed speeds do arise from 
the plasma bulk motion, then the above measures place strong con- 
straints on its velocity and the angle of the jet to the line of sight. 
Biretta et al. (1995) found that (j> < 43°. Assuming that the ap- 
parent speeds are ~ 6c then <f> < 19° (Biretta et al. 1999) (further 
allowed combinations of F and 4> are shown in their Table 3). As 
a result we use F = 10.0 ± 5.0 as this encompasses a large pro- 
portion of the range of observed superluminal speeds, along with a 
corresponding <f> — 15 ± 5°. 

Dodson et al. (2006) measure a spectral index of —1.1 ± 0.4 
between 1.6 and 4.8 GHz, and -0.6±0.4 between 4.8 and 15 GHz. 
Their calculations of the spectral turnover gives ~ 20 GHz. The 
core flux densities used here have been measured at 5 and 15 GHz. 
These are close to the spectral turnover frequency as estimated by 
Dodson et al. (2006) and are thus presumably in the self-absorbed 
region of the spectrum, where Equation 3 is applicable. 



3.2 Perseus 

Studies of the core of 3C 84 have also resulted in a range of mea- 
sured apparent motions. Kellermann et al. (2004) using 15 GHz 
VLB A observations find f3 = 0.2 ± 0.1. Walker et al. (1994) find 
that the jet speed and orientation angle are not tightly constrained, 
P = 0.3 - 0.5 and = 30° - 55°. As in M87, some compo- 
nents appear to be much slower, (3 ~ 0.05 (Dhawan et al. 1998), 
and others are faster, /? = 0.7 ± 0.1 (Marr et al. 1989). Krichbaum 
et al. (1992) estimate an angle to the line of sight of ~ 10° and a 
corresponding j3 > 0.68. We use the values presented in Walker 
et al. (1994) as these fall towards the middle of the observed range 
09 = 0.3 - 0.5 corresponds to T = 1.10 ± 0.05). 

The core of 3C84 in Perseus has a complicated spectral 
index structure with an inverted spectrum at lower frequencies 
(ai. 5-5.0 GHz ~ 1, Unwin et al. 1982; O'Dea et al. 1984), flat- 
tening towards higher frequencies ccio-iooghz ~ (O'Dea et al. 
1984). The spectra in Readhead et al. (1983) also show a decrease 
of flux density toward lower frequencies, though there are no mea- 
surements at less than 10 GHz. The spectral index varies across 
the core and so an exact value is difficult to obtain (e.g. Unwin 
et al. 1982; Readhead et al. 1983; Vermeulen et al. 1994; Krich- 
baum et al. 1992). Silver et al. (1998) indicate that the spectrum 
may indeed be flat in the range 1 — 10 GHz (see their Fig. 5). Our 
core flux for Perseus is at 15 GHz and so it is reasonable that Equa- 
tion 3 is applicable. 
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Table 2. : Core fluxes and jet properties 



Name 


Date 


Flux Density (Jy) 


Size (mas); 6 a , 9b 


pa 


!/(GHz) 


Reference 


M87 


Reynolds' 3 


0972/0373 


1.00 ±0.10 


0.70 ±0.07 


0.70 ± 0.07 


3±1 


5.0 


1,2,3,4 


Kovalev-Score 


050203 


1.39 ± 0.07 


0.41 ±0.04 


0.27 ±0.03 


10 ±5 


15.0 


5,6 


Kovalev-S? T 

Unres 


050203 


0.73 ± 0.04 


0.41 ±0.04 


0.27 ±0.03 


10 ±5 


15.0 


5,6 


Lister05-7c 


050203 


0.98 ± 0.05 


0.37 ±0.02 


0.14 ±0.01 


10 ±5 


15.0 


6,7 


Scott-S c 


201297 


0.18 ± 0.03 


0.20 ±0.02 


0.20 ±0.02 


10 ±5 


5.0 


6,8 


Perseus 


Kovalev-Score 


010303 


3.63 ± 0.18 


0.66 ±0.07 


0.21 ±0.02 


1.10 ±0.05 


15.0 


5,9 



NOTES: All fluxes and sizes are model fits to the emission. The beam for Lister05 is 1.23 X 0.55 mas. Date is in DDMMYY or MMYY/MMYY format. 
a The r is for the bulk motion of the jet. b a = 0.5 for this core flux. c The angular sizes are an upper limit for these measurements. REFERENCES: 1. 
Pauliny-Toth et al. (1981); 2. Biretta (1993); 3. Biretta et al. (1995); 4. Reynolds et al. (1996); 5. Kovalev et al. (2005); 6. Biretta et al. (1999) 7. Lister & 
Homan (2005); 8. Scott et al. (2004); 9. Walker et al. (1994). 
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Figure 1, Constraints on the B — n plane from Equations 2 and 4. The 
expected values for the number density for light and heavy jets from the 
estimated Lk are shown and correspond to the circle and square respec- 
tively. The "intersection point" (triangle) is the number density required 
from radio observations. This schematic corresponds to Kovalev-Sfjnres 
from M87, and for clarity the uncertainties are not shown. 



4 RESULTS 

Fig. 1 shows, for a representative case, the constraints on the mat- 
ter content of jets placed by Equations 2 and 4 in the B — n plane. 
The lines labelled "Eqn. 2" and "Eqn. 4" exclude regions of the 
B — n plane. The "intersection point" (triangle) is thus the mini- 
mum number density possible in the jet from observations, being 
the combination of the two above equations. The expected number 
densities inferred from the estimated Lk for pure electron-positron 
(the circle and line labelled "n cc ") and electron-proton (the square 
and line labelled "n cp ") jets are shown. In this representative case 
an electron-proton jet is excluded, and hence the jet is electron- 
positron. Note that if the estimate on the magnetic field were to 
increase then the intersection point would move to the right, even- 
tually allowing the possibility that the jet is electron-proton. When 




Figure 2. The relative number densities in the jet Circle for electron- 
positron and Square for electron-proton. Triangle makes the intersection 
of the two limits in the B — n plane and represents a lower limit. 



this occurs, no constraints can be placed on the matter content of 
the jets. 

The results are shown in Fig. 2, for both M87 and Perseus. 
In order to compute Lk we average the energies and timescales 
of the bubbles within each cluster, as the differences in the results 
between the different bubbles are very small. 

Fig. 2 reports the intersection points corresponding to the 
n — B plane for each useful set of data. In order to check the ro- 
bustness of the conclusions, we also explicitly estimated errors on 
the inferred values. In particular, the reported uncertainties are es- 
timated using a Monte-Carlo algorithm, assuming the input uncer- 
tainties have a Gaussian distribution which take into account range 
of values measured for each parameter. The source of uncertainties 
which have been accounted for are those parameters which have 
an uncertainty quoted in Tables 1 and 2. If no uncertainty has been 
given in the literature, we have adopted a conservative value of 10 
per cent (a is assumed to have no error, but the effect of different 
values is discussed in Section 6.2). The error-bars shown on Fig. 2 
are 1 — a errors derived from this algorithm. These uncertainties 
are only applicable to this model for the core radio emission. 

Although the different datasets imply different values of den- 
sities, globally they provide a similar result; that the relativistic jets 
in these two sources are likely to be electron-positron. This analy- 
sis indicates that for M87, including the core flux used in Reynolds 
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et al. (1996), the jet is electron-positron. There is, however, only 
one measurement for the Perseus cluster (from Kovalev et al. 2005) 
which places clear limits on the matter content of the jet in 3C 84. 

As a result of the large uncertainties in the measurements of 
the jet angles and bulk velocities and the disparity between the val- 
ues for M87 and Perseus we check if the conclusions change sub- 
stantially using different V and (f>. We recalculate the results for 
M87 with T = 1.1 ± 0.05 and <J> = 35 ± 5° and for Perseus with 
r= 10 ± 5 and 0= 10 ±5°. 

Relative to the number densities required by the estimated Lk 
(which themselves rise by a factor xlO 2 — 10 3 ) the intersection 
points (red triangles) for M87 fall. The results still favour light jets, 
though are less definite. In Perseus, as the jet angle and velocity 
were altered in the other sense, the number densities corresponding 
to Lk decrease and the limits become more consistent with light 
jets. 

Thus, even though we have chosen values for T and cj> which 
are drastically different from the average of the measured values 
the results still indicate that both jets are light, showing that the 
conclusions do not depend critically on these values. 

Other clusters were also analysed, including A2199 (3C 338) 
and Hydra A, but the radio observations of the core did not have 
sufficient resolution (the model is very dependent on the angular 
size of the core). The resultant magnetic field estimates were there- 
fore high, giving very low estimates on the number densities in the 
jet. Therefore no useful constraints could be placed on the matter 
content. 
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Figure 3. n — 7 m m parameter space. The thick lines report the constraints 
imposed by the four conditions: lower limits from SSA (solid black line); 
upper limits from SSC (dashed black line); jet kinetic power from electron- 
positrons (green line) and electron-proton (blue line); upper limit from the 
Larmor frequency v B (vertical magenta line). The thin lines define the up- 
per and lower uncertainties on each constraint as derived via the Monte 
Carlo simulations. The figure illustrates the case for the Kovalev-Score 
dataset for M87. 



5 7min LIMITS 

We have so far assumed that the low-energy cut-off in the energy 
distribution of leptons occurs at 7 m i n ~ 1. Celotti & Fabian (1993), 
Ghisellini et al. (1992) and Reynolds et al. (1996), however, find 
that a low-energy cut-off at 50 MeV (7 m i n ~ 100) in the lepton 
energy spectrum is required for both heavy and light jets, consistent 
with polarisation measurements. 

As 7 m in rises then the number densities will fall as fewer par- 
ticles will carry the same energy (the average energy per particle 
will rise). In the assumption that the protons remain cold, one can 
even reach the situation where the average lepton energy is close to 
the proton energy (i.e. ( 7 ) ~ — ). It is clear that when approach- 
ing this value, the expected lepton number densities for the two jet 
types differ by a smaller and smaller factor and so, given the large 
uncertainties of the current measurements, by increasing 7 m i n the 
possibility of drawing any conclusions on the matter content be- 
comes increasingly small. 

5.1 Synchrotron and Synchrotron Self-Compton constraints 

Upper limits on 7mm in the core, however, can be set. 

A first constraint on maximum values of 7mm is set by the 
very same self-absorption model adopted above, as, for it to be self- 
consistent, the lepton distribution should extend at least to energies 
of the leptons responsible for the emission at v m . In other words, 

Further upper limits can be also obtained from synchrotron 
Self-Compton (SSC) constraints. Namely, the X-ray emission ex- 
pected from the self-absorbing region via SSC, Sx should not ex- 
ceed the observed X-ray flux, where 

Sx= 2aS m na Tn ^ n M° h M (6) 
t{a) V^x/ V^m/ 



(e.g. Reynolds et al. 1996; Marscher 1987). Here <tt is the Thom- 
son scattering cross-section, and t(a) is tabulated in e.g. Ghisellini 
et al. (1992). The break frequency/upper cut-off in the synchrotron 
spectrum, Ub, is taken to lie at ~ 10 lj Hz (Meisenheimer et al. 
1996) and fx is the X-ray frequency. It is of course possible that 
the X-ray emission is not entirely the result of the SSC process. 
Some may be synchrotron emission (e.g. Marshall et al. 2002), in 
which case S% is an upper limit to the SSC flux. 

These limits are, as for the SSA, dependent on wy^ n , while 
the kinetic luminosity depends on 71(7) . It is thus necessary to com- 
pare the above constraints in the n — 7mm parameter space for 
7min > 1. Figure 3 reports the results for a representative case 
(referring to the dataset Kovalev-Scorc of M87), each thick line re- 
ferring to the upper limit from SSC (dashed black), the lower limit 
from SSA (solid black), the values corresponding to the inferred 
Lk for pair (green) and electron-proton (blue) jets, and the upper 
limit on 7 m i n from vb (magenta). 

Taken at face value, the implication is that all the conditions 
can be satisfied by either a pair jet with 1 < 7 m in < 100 or a heavy 
jet for the limited range 50 < 7 m in < 100. Qualitatively similar 
inferences can be derived for M87 in the cases of Kovalev-Sunres 
and Reynolds (fully in agreement with the conclusion reported by 
Reynolds et al. (1996)), while the Lister- I c is just marginally con- 
sistent with pair jets with 10 < 7mm < 150 and Scott-5 1 does 
not formally allow a consistent solution. These last two core flux 
measurements are those in which the minimum number density 
fell above that allowed for an electron-positron jet in Fig. 2. For 
all cases the SSC upper limits do not give interesting constraints. 
On the contrary the Perseus X-ray data (SSC) strongly constrain 
the possible solutions to light jets with 7 m in < 10 or heavy jets 
with 10 < 7min < 50. Overall, 7 m in < 10 would imply electron- 
positron jets. For higher values, typically 10 < 7 m in < 100 both 
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heavy and light jets are allowed by the data. These 7 m i n limits are, 
on average, not inconsistent with those put forward by Fabian et al. 
(2002). It should be nevertheless stressed the (in some cases im- 
plausibly) tight range of values of 7 m ; n were the jets heavy. 

The errors on the source parameters, however, limit the pos- 
sibility of drawing any firm conclusion at this time. This can be 
seen in Fig. 3, where the errors on the above constraints, again de- 
rived via the Monte Carlo simulations, are reported as thin lines. 
It should also be noticed that the condition imposed by could 
be in principle more relaxed if the upper limit on B has been over- 
estimated because of the limited resolution of the observations, as 
7min oc B^ 1 ^ 2 oc 6~ 2 (Equation 2). 

We conclude that currently this method cannot definitely rule 
out electron-proton jets in favour of electron-positron jets, although 
the former case would require a rather narrow range of values of 
possible 7mi n . Better quality data (angular resolution, jet speed, jet 
direction, contribution of SSC emission in the X-ray band) and/or 
other constraints on 7 m i n could impose firmer conclusions. 



6 DISCUSSION 

In the following we will summarise and discuss the caveats of the 
above estimates, from both observational and modelling points of 
view. 

At a general level, the conditions on small scales in the jet 
may not correspond to those on large scales, from conservation of 
particles, momentum and energy alone. There may be entrainment 
of ambient material by the jet, dissipation, turbulence or Poynting 
flux energy transport. 

The jet of M87 has filamentary structure on all scales where 
it is well imaged. There is no obvious distance from the nucleus 
where the jet disappears. If there are any inhomogeneities on small 
scales in the jet, then the synchrotron self-absorption estimates 
could be affected. 



6.1 Variability & Kinetic Luminosity 

These calculations assume that the jet is of constant power, corre- 
sponding to that observed now. This is not necessarily true. As the 
synchrotron flux of the core has been used to determine the matter 
content of the jet, any flux variability would change this (both den- 
sity and magnetic field), particularly so if the size of the emitting 
region appears the same (see Equation 2). It would, however, have 
no immediate effect on the bubble, and hence the inferred kinetic 
luminosity would remain the same. 

As B oc 5 m 2 and nB 3/2+a = nB 2 ' 2 = const (Equation 4), 
then n oc S^ 4 (for a constant 6 and a = 0.7). So a small change in 
the core flux changes the minimum possible number density in the 
jet, even by orders of magnitude. 

To see by how much the core would have to vary for the 
number densities in M87 to be consistent with electron-proton jets 
the core flux for the Kovalev-SWres and Scott-S* values were de- 
creased until the intersection point is consistent with the electron- 
proton values. No other values were changed except the flux. The 
Kovalev-Sunrcs core flux only has to change by a factor of 4 for the 
number densities to be consistent with electron-proton jets. In the 
case of Scott-5 1 the variability has to be around a factor of 8. 

Indeed, the flux of 3C84 steadily increased between 1960 and 
1989 (Marr et al. 1989; O'Dea et al. 1984), by up to a factor of 
4. More recent observations in the optical and X-ray show that the 



nuclear flux varied by less than a factor of 3 between 2000 and 2003 
(Perlman et al. 2003). 

This is a major problem with the method utilised in this work. 
In principle, there is no temporal link between the current kinetic 
luminosity of the jet and that calculated from the bubble. The ra- 
dio core and the jet fluxes may vary over the course of decades 
whereas the bubbles have timescales of ~ 10 7 yr. However re- 
cent work by Allen et al. (2006) found a correlation between the 
Bondi accretion power and the kinetic luminosity of jets in nine 
nearby elliptical galaxies. The bubble powers are measured over 
a timescale of ~ 10 7 yr and the inflow through the Bondi radius 
takes ~ 10 4 — 10 5 yr. Therefore, the resulting correlation suggests 
that these powers might be comparatively steady over ~ 10 7 yr. 

Further to the discussion on the bubble powers (Section 2.3) 
is the issue of any energy in the form of cosmic rays and magnetic 
fields present at the centres of clusters. Sanders et al. (2005) show 
that the pressure of these components could contribute at least 30 
per cent in the region of the Perseus bubbles. As a result the energy 
content of the radio bubbles may be larger than just what thermo- 
dynamic calculations imply, requiring a larger number density of 
particles in the jet. There may also be losses as the jet travels out 
from the central AGN into the lobe. The kinetic luminosity inferred 
from the bubble expansion would be an under estimate of the initial 
kinetic luminosity of the jet. As such the expected number densities 
may be lower limits. These effects could cause core fluxes which 
currently exclude heavy jets to become consistent with them, but 
only if these effects are (unexpectedly) large (x ~ 100). 

6.2 Parameter Effects 

As was presented in Sections 3.1 and 3.2 there is no clear, exact 
consensus on the value for all the observational parameters which 
are used in this model. The true jet velocity and angle to the line 
of sight are still uncertain, and the range in the currently inferred 
values is large. 

The jet velocity affects the expected number densities (Equa- 
tion 5), Lk oc nT 2 and the magnetic field estimate (Equation 2), 
B oc 8. Therefore, if the true jet velocity were faster, then the ex- 
pected number densities required to create the observed bubbles 
would fall. If so, the two core flux measurements for the Perseus 
jet not presented in the main results (from Kovalev et al. 2005 and 
Scott et al. 2004) could also imply electron-positron jets (see also 
Section 4). 

The angle to the line of sight changes the magnetic field cal- 
culated in the jet, via the Doppler factor 8. The change (increase) 
in S becomes greater as the angle to the line of sight (f> decreases. 
As a result the upper limit on the magnetic field increases and the 
"intersection point" decreases, making the constraints on the mat- 
ter content less restrictive. We have, however, taken into account 
the large uncertainties on F and cf> in our calculations of the errors 
associated with the estimated number densities. 

The spectral index of the optically thin spectrum, a, has in- 
stead very little effect on the results in the range 0.5 < a < 1.0. 
As such the conclusions remain unchanged over a large range in 
a which is a parameter whose value we have assumed rather than 
measured. 



6.3 Dimension of the Self- Absorbed Core 

The angular size of the core has been taken as the diameter of the jet 
at the point when the jet becomes self-absorbed; however the core 
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may not be resolved and so this may be an upper limit. Although 
we have used the most recent high resolution data, the actual emis- 
sion region might be smaller than the size of the radio telescope 
beam, therefore we would have over-estimated the angular size of 
the source. As B oc 9* (Equation 2), the effect of any difference 
between the true and observed value would be greatly amplified in 
the estimates of B and n. Indeed, the sources with large estimates 
of magnetic field (e.g. A2199, Hydra A) were also those furthest 
away, implying the resolution may not have been high enough to 
fully resolve the emitting regions. 

Another caveat is that the angular size of the core may 
also represent the entire projected length of the self-absorbed jet 
(Marscher 1987). However, the most recent observations of M87 
have a resolution of 0.1 milliarcsecond (Krichbaum et al. 2004), 
which corresponds to tens of light days. This means that at least for 
the closest objects the core radii are unlikely to be blends of knots 
of emission, masking the true core emission. 

6.4 Comparison with other studies 

Ghisellini et al. (1992) used spectral information and VLBI maps 
from 105 sources along with an SSC model to obtain estimates on 
the bulk Doppler factor and emitting particle number density. Up- 
per limits on the pair number density were placed from the num- 
ber expected to be able to survive annihilation from the central re- 
gion where they might be produced, and lower limits from the ob- 
served synchrotron emission. They found that a low-energy cut-off 
of 7min ~ 100 in the lepton energy distribution was required for 
both light and heavy jets. 

Using the same sample, Celotti & Fabian (1993) found that 
the estimated kinetic power on parsec scales was similar to that 
inferred from radio lobes, concluding that energy can be efficiently 
transported over several orders of magnitude in length scale along 
the jet. Combining limits on the total kinetic power and particle 
number density flux, their results favoured heavy jets, also with a 
low-energy cut-off at around 7 m i n w 100. 

From X-ray observations of blazars associated with optically 
violent variable quasars, Sikora & Madejski (2000) place con- 
straints on the pair content of the radio-loud quasar jets. They ex- 
clude both pure light jets, as these over-predict the soft X-ray flux, 
and pure heavy jets as these predict too weak non-thermal X-ray 
emission. Although the pair number density is larger than the pro- 
ton number density, the jets would be dynamically dominated by 
the protons. 

Circular and linear polarisation observations of jets can be 
used to constrain the low-energy particle distribution, the magnetic 
field strength and the particle content. Depending on whether this is 
intrinsic to the synchrotron emission or produced by Faraday con- 
version of linear polarisation to circular, different limits can be set 
on the low energy particle distribution (e.g. Homan 2005). 

Wardle et al. (1998) measured circular polarisation from 
3C 279 and, considering that most likely this results from Fara- 
day conversion, set an upper limit 7 m i n < 20, which would be 
evidence for electron-positron jets in this source. Further observa- 
tions of PKS 0528+134, 3C 273 and even 3C 84, detected circular 
polarisation, but no limits on 7 m i n have been determined (Homan 
& Wardle 1999, 2004). In any case, more recently Ruszkowski & 
Begelman (2002) showed that the observations of 3C 279 could be 
consistent with both types of jet, by arguing that the linear and cir- 
cular polarisations observed could, depending on the field configu- 
ration, be consistent with different plasma compositions. Therefore, 
no strong conclusion can (yet) be drawn on such measurements. 



While we have no robust explanation on the discrepancy be- 
tween results supporting electron-proton dominated jets and our 
findings, one intriguing possibility is that jets associated with pow- 
erful sources are energetically dominated by a proton component, 
while low power radio sources, such as M87 and Perseus, are pre- 
dominantly composed of an electron-positron pair plasma. 

6.5 Magnetic Field Energy 

In general, it is possible that jets are magnetically dominated, with 
negligible contribution to the energetics from particles (e.g. Bland- 
ford 2002). Sikora et al. (2005) argued instead that conversion from 
Poynting-flux dominated to matter dominated jets takes place very 
close to the black hole. 

From arguments of energy equipartition between the emitting 
leptons and the magnetic field, B oq ~ 10 -6 to 10~ 3 G are es- 
timated depending on the jet scale. Studying M87, Stawarz et al. 
(2005) find a lower limit B > 3 x 10~ 4 G in the brightest knot 
of the jet at ~ 2 kpc from the core, concluding that it is likely that 
this represents a departure from equipartition, such that the knot is 
at least marginally magnetically dominated. 

Any Poynting flux would not be easily detectable from the ra- 
diative properties of the jet (except possibly from a population of 
particles with a narrow energy distribution, as expected following 
acceleration from reconnection events) and because of that no con- 
clusion could be drawn from the current data. 

In this work we have assumed that the magnetic field does not 
dominate the energy content. 



7 PARTICLE CONTENT OF RADIO LOBES 

The results imply that the jets in M87 and 3C 84 are light (electron- 
positron), at least at their base. Fabian et al. (2002) quantify the 
particle content of the radio bubbles in the following way. By com- 
paring the pressure inside and outside of the bubble, any extra par- 
ticles beyond those inferred from the radio emission can be deter- 
mined. This method gives estimates on k/f, where / is the vol- 
ume filling factor of the relativistic plasma and k accounts for the 
extra particles. Dunn & Fabian (2004); Dunn et al. (2005) find 
1 £ k/f < 1000 (for M87 k/f ~ 10 and k/f ~ 500 - 1000 
for Perseus) and that the radio lobes in these clusters are not in 
equipartition but are particle dominated. 

If the jets start as a pure electron-positron plasma, then for pro- 
tons to be present in the lobes, some entrainment of material has to 
occur. This is likely to be a stochastic process and highly dependent 
on the environment surrounding the jet within the inner few kpc. 
Powerful sources are unlikely to pick up much material as their jets 
would clear out the regions surrounding them and weaker sources 
would be expected to pick up more material. Entrainment may be 
the key reason why the latter sources appear as FR I, as indeed de- 
tailed modelling indicates (e.g. Laing & Bridle 2002). In principle, 
entrainment could slow down a jet to sub-relativistic speeds. The 
amount of entrainment, however, especially on large scales where 
the jets are no longer fully collimated, is uncertain. The exact shape 
of the low-energy lepton spectrum will also affect k/f, as it will the 
matter content of the jet. 

De Young (2006) suggest that the fact that k » 1 (for / = 1) 
in some clusters results from a population of cold protons present 
in the jet. This is different to the jet matter content suggested here. 
This solution does have problems as the accompanying relativistic 
electrons may cause jet to decollimate unless magnetic confinement 



© 0000 RAS, MNRAS 000, 000-000 



8 Dunn, Fabian & Celotti 



occurs. Another solution De Young (2006) suggest are Poynting- 
flux dominated jets, though these also have associated difficulties, 
including whether they can produce the FR I morphologies ob- 
served (e.g. Komissarov 1999; Leismann et al. 2005). 

The radio lobes of M87, Cygnus A and Fornax A appear very 
filamentary and so the filling factor of the radio emitting plasma 
may vary, which may have bearing on the matter content of the 
jet. The net effect is not obvious as it would depend on the actual 
structure of the relativistic component and its possible confinement. 
While clearly the kinetic luminosity and synchrotron emission in 
the optically thin regime only depends on the total particle number, 
the synchrotron self-absorption estimates are also affected by the 
density inhomogeneities along the line of sight. 



8 CONCLUSIONS & OUTLOOK 

If the conditions observed at the base of the jet are typical over the 
timescale for bubble creation, and all the caveats discussed in the 
previous section hold, then the jets in M87 and Perseus are likely to 
be dominated by an electron-positron plasma. Other radio sources 
in clusters analysed are too distant to afford the resolution neces- 
sary to measure the synchrotron flux from the base of the jet. 

Despite the significant potential offered by estimates of the jet 
power from its interaction with the ICM, this method is currently 
limited to determining the matter content of extra-galactic radio 
jets in sources in sufficiently nearby clusters of galaxies. A con- 
certed effort in monitoring the cores of nearby AGN at a variety of 
frequencies will help in constraining the jet matter content. Clear 
identification of the cores, measurements of the jet velocity and 
limits on the variability will aid in improving the parameters and 
assumptions involved in these calculations. Future perspectives in- 
clude high accuracy measurement of resolved radio cores via the 
space VLBI, and higher frequencies observations (though the flux 
should not be dominated by optically thin emission). 

A definitive result on the matter content of jets remains, how- 
ever, elusive for now. 
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